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ABSTRACT 

Two hundred and eighty-five hours of solar data obtained from 
the University of Iowa 2 cm radiometer during 1968 — 1969 were analyzed 
for evidence of the 5 min periodic chromospheric oscillations. A 
power spectral analysis of the data failed to show any statistically 
significant (> 9 confidence) periodic activity in the frequency 
range 1 — 15 mHz. No correlation was found between 2 cm periodicities 
and other' forms of solar activity at several different chromospheric 
heights in H ,■ soft solar x rays (2 — 12 A), and several microwave 
frequencies (3 — 15 GHz). 

A small shift in power from low to higher frequencies in the 
power spectrum of the 2 cm data was found to be correlated with H 
and x-ray activity. This power shift is attributed to a relative 
increase in chromospheric turbulence at altitudes common to H , x-ray 
and 2 cm emission. 

A statistical analysis of previous works reporting evidence for 
the oscillations at microwave frequencies indicates that confidence 
in these previous results is marginal. 

A model for chromospheric oscillation bursts in quiescent 
supergranules deduced from Bhattacharyya (1972) is presented. The 
model is an exponentially damped wave train burst with a 15 min mean 
damping time and a mean burst period of 50 min. Using the statistical 
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properties of power spectra, an expression for the minimum detectable 
signal-to-noise ratio is derived and applied to the model using the 
most significant results of experiments to date to yield: 

(1) an estimate of ~ 1.5 X 10 as the upper limit for the ratio of 
the oscillation mean peak-to-peak to ambient temperature at 
chromospheric heights, 

(2) an estimate of ~ 0.02 as the upper limit for the mean fraction 
of the quiescent solar surface experiencing oscillations, and 

(3) estimates for the minimum amount of data required for confident 
detection of the periodicities. In an ideal experiment utiliz- 
ing the NRAO Kitt Peak 3 6 ft antenna at a wavelength of 1 mm 
and assuming the present model for oscillation bursts and rea- 
sonable equipment and data parameters, ~ 1.25 hr of data are 
calculated to be the least amount required. 
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I. INTRODUCTION 

The existence of characteristic oscillations in the solar 
atmosphere has been known for some time. Following the first theo- 
retical treatment of the problem by Whitney (1958) the exper im ental 
verification of the 5 min period oscillations was reported by 
Leighton et al. (1962) using a Doppler spectroheliograph technique. 

A short time later Evans and Mitchard (1962) and Howard (1962) 
extended the observational techniques of detecting the periodicities 
by using, respectively, time series analysis of high dispersion 
spectra and photoelectric detection. Since that time many investi- 
gators have succeeded in extending the range of optical observations 
using a variety of techniques over a wide range of spectral frequen- 
cies. Noyes (1967) gives a comprehensive summary of early observa- 
tions and the discovery of the solar supergranulation and its rela- 
tion to the chromospheric network, and reviews the geometrical model 
which resulted from those observations. 

An important finding was that the oscillations existed at 
heights in the solar atmosphere extending into the lower chromosphere. 
Generally speaking a peak in the velocity spectrum occurred at a per- 
iod of about 300 sec for spectral lines corresponding to low photo- 
spheric altitudes and moved toward shorter periods as observations were 
made at higher altitudes (Noyes and Ieighton , 1963)* Th e intensity 
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spectrum followed the. velocity spectrum closely, though it was more 
pronounced at high photospheric altitudes where the atmosphere is 
nearly adiabatic than at low photospheric altitudes where it is 
nearly isothermal (Noyes and Leighton, 1963). The observation of the 
oscillations in the H velocity spectrum (Elliot, 1969 ) and in the 
1^, intensity spectrum (Bhatnagar and Tanaka, 1972) showed that the 

oscillations persist to heights characteristic of H line formation, 

3 

on the order of 5 * 10 km. Since this is in the neighborhood of 
the region where millimetric and centimetric radiation are observed 
to originate it is reasonable to ask whether or not the observed 
fluxes at these radio frequencies are modulated by the pressure 
fluctuations associated with the oscillations. 

Reported radio observations of the oscillations have been far 
fewer in number than optical observations. Yudin (1968) reported 
having observed the 5 min period intensity modulation of 3-3 cm 
radiation. This same period, along with several other components, 
was also reported by Durasova at al. (1971) at the same frequency 
and apparently by using the same instrument. Simon and Shimabukuro 
(1971) found an oscillatory component at 5.6 mHz corresponding to 
a period of 180 sec, in their study of 3.3 and 3-5 mm solar data. 
Longer period components of from 30 to 60 min were investigated by 
Kobrin and Korshunov (1972) at a wavelength of 3.1 cm and by 
Kauffman (1972) at 4.3 cm. Gotwols (197l) showed that short period 
oscillations can exist after flares in the decimetric wavelength 


region. 
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Establishing the degree to which the oscillatory phenomenon 
exists at chromospheric and higher altitudes is important for a 
number of reasons, perhaps chief among which is determining how large 
a role it plays in heating the solar corona. A search was instituted 
for periodicities in the range 1 — 15 mHz using data taken in 
1968 — 1969 on the University of Iowa 2 cm radiometer. The results 
of the search and a review of the findings of others, ^with particular 

I 

emphasis given to the statistical aspects of the resulting power 
spectra, are the subject of this thesis. 

Section II describes the equipment used to gather the data and 
the data reduction and digitization procedures. The data activity 
classification scheme and power spectrum analysis are described in 
section III. Section IV is devoted to a discussion of the results 
of the power spectrum analysis. The absence of detectable periodi- 
cities is compared with the findings of others and an attempt is 
made to assess the statistical reliability of the respective results. 
Several statistical parameters related to the power spectrum detect- 
ability of a signal embedded in noise are derived in section V. 

These parameters are applied to an oscillation burst model to yield 
upper limit estimates for oscillation amplitudes and the mean frac- 
tion of surface area experiencing oscillations in the quiescent solar 


atmosphere. 
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II. EQUIPMENT DESCRIPTION AND DATA HANDLING 

i 

Data used for this study were recorded from the 2 cm radio- 
meter operated by the University of Iowa at the North Liberty Radio 
Observatory (NLRO) (Wende, 1969 ). The system has a 4 ft diameter, 
equatorially -mounted aluminum paraboloid equipped with a dual linear 
polarized Cassegrain feed. The antenna half-power beamwidth is 1° so 
the entire sun is viewed. The receiver operates at a center frequency 
of 15.375 GHz with an IF bandwidth of 20 MHz. 

The radiometer is of the Ryle-Vonberg type. Antenna noise 
is balanced against the noise from, a gas discharge tube-precision 
attenuator, resulting in a stability dependent only upon the gas 
discharge tube stability. The servo motor is driven at 100 pps, 
yielding a 10 msec response time for a 4 K change in antenna tempera- 
ture. A potentiometer is connected to the servo motor shaft, pro- 
viding an output voltage according to the position of the attenuator. 
This output voltage is used to drive a voltage controlled oscillator 
(VCO). For a more complete description of the system see Wende 
( 1968 , 1969 ) and Collins ( 1967 ). 

During operation the radiometer was calibrated daily using a 
procedure which entailed inserting a known servo error voltage from 
a precision potentiometer. From the manufacturer's calibration 
curves one could then determine the calibration curve of the, antenna 
temperature T^ versus the VCO frequency as shown in Figure 2. 
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The VCO output was recorded on magnetic tape along with a WWV 
calibrated, 1 kHz NASA 36-bit time code timing reference. The record- 
ed timing reference was later used to trigger the A/D digitizing 
equipment and provided an absolute standard against which to measure 
the integrated VCO frequency. It was thus insured that any effects 
due to wow and flutter present in either the recording or playback 
systems were effectively compensated for. 

A/D conversion was carried out by separating the magnetic 
source tape output into its time code and data components and using 
the time code reference to control the data sampling system, as repre- 
sented schematically in Figure 3. The sampling system had a maximum 
time resolution of 100 msec and basically consisted of a frequency 
counter with a time base set equal to the desired sample period, which 
in this case was 5 sec. The output from this system therefore repre- 
sented a 5 sec integrated average of the radiometer VCO frequency. 
Integrated samples were recorded sequentially on paper tape. One 
time sync was also punched every 10 min, containing the start time 
of the 10 min period, sampling period, and supplemental information 
(such as weather information for that data day) that could be coded 
into the system by means of switches on the front of the sampling 
system. The digitized, integrated VCO frequency samples were con- 
verted to antenna temperature by means of the calibration curve in 
Figure 2 and stored on magnetic tape for later computer analysis. 
Extensive checks were made during the conversion process to detect 
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spurious data points which were replaced by values interpolated from 
good points on either side. The fraction of data points that had to 
be corrected in this way was less than 0.01 for the entire set of 
digitized data and is estimated to be less than 0.002 for the data 
selected for power spectrum analysis. 

Estimates can be made for the receiver noise (Krauss, 1966) 
and the effective temperature resolution of the system. 

For the receiver noise 



T s /B At 


where 

AT = rms receiver noise , 
r 

T g = system temperature 

= T A + [(10 log 1Q NF) - 1] T q , 

T^ = antenna temperature , 

NF = system noise figure w 11 Db, 

T q = ambient temperature *=* 500 K, 

B = receiver bandwidth = 20 MHz, and 
At = sample integration time. 

For a typical day T^ ~ 1000 K. In order to achieve a reason- 
able compromise between statistical reliability and frequency resolu- 
tion of the computed power spectra At was chosen to be 50 sec. The 
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5 sec integrated T. averages were therefore averaged in groups of 

.H. 

six to yield an effective data sampling period and integration time 

of 30 sec. These parameters then give a system noise of AT « 0.l6 K. 

r 

Although the stepping resolution of the system is approximately 
4 K, for a noisy signal with an rms deviation appreciably larger than 
this the effective temperature resolution of the system for a sample 
period of AT is given by 


AT 


AT 


eff 


step_ 


v/f 


step 


AT 


where 


AT 


eff 


AT 


step 


step 

At 


= effective stepping resolution, K; 
= system stepping resolution, K; 

= stepping frequency = 100 Hz; and 
= data sample period = 30 sec. 


Hence AT ^ 0.07 K. 
el 1 

There is an additional source of noise due to the uncertainty 
of the sampling system frequency counter of ± 1 count per data sample. 
For a VCO frequency of 600 Hz (T a w 1000 K) this results in an 
uncertainty of 1/3000, and for six averaged 5 sec samples corresponds 
to a probable error factor of (l/3000) X (l//£) = 1.5 X 10~\ From 
Figure 2, this factor corresponds to AT^ « 2 K. 

Thus the largest source of artificially-generated noise in 
the reduced antenna temperature data is due to that introduced by 
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the digitization process. However the rms noise of the quiescent 
sun about the mean at any given time is observed to be on the order 
of 10 — 20 K for a data sample integration time of 30 sec, or nearly 
an order of magnitude larger than AT^. We therefore conclude that 
all but a small fraction of the noise present in a plot of versus 
time is true solar centimetric noise. 
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III. DATA ANALYSIS 


Of the more than 200 days of data taken during 1968 — 1969* ^0 
days, representing those of the best quality, were selected for study. 
The criteria for selection were that the data must be free of dis- 
continuities and irregularities, weather conditions must be clear, 
and a minimum of 4 h of continuous data must exist in a given record. 
For a given day records as long as possible were chosen to be free 
of any large flares that might force the data to depart signifi- 
cantly from stationarity, though no attempt was made to choose data 
completely free of all activity. The total amount of data repre- 
sented in all records was greater than 285 h, making the average 
record length about 7 h. The 5 sec data samples were averaged in 
groups of six to yield samples corresponding to 50 sec integrated 
antenna temperature averages. A typical day of data with 50 sec 
time resolution is shown in Figure 4 . 

The standard Blackman-Tukey (1958) algorithm with Hamming 
smoothing was used to compute the power spectrum for each data record. 
The maximum lag M in the autocorrelation function was set at 25 min 
for each record so that all of the individual spectra have the same 
frequency resolution Af = l/2M. Because of the large amount of 
power present at low frequencies the data was pre-whitened prior to 
computation to prevent the low-frequency power from leaking to higher 
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frequencies via the side lobes in the spectral window. This was 
accomplished by means of a digital high- pass filter constructed as 
follows . 

For n successive passes through the data with the elementary 
low-pass filter function 

V = 1 x i + 5 < x i-i + x i + i> ’ 

Holloway (1958) showed that the power transfer function is 
R(f) = cos 2 n (TTfAT) 

If the resulting smoothed data is subtracted from the original 
unsmoothed data set the difference represents the original data 
filtered by a high-pass filter with a power transfer function equal 
to the complement of R(f), 

R / (f ) = 1 - R(f) 

The antenna temperature data were filtered using this technique with 
n = 10. The resulting transfer functions R(f) and R / (f) are plotted 
in Figure 5 • After pre-whitening was completed the data were nor- 
malized to zero mean and unit standard deviation and the power 
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spectrum computed. Post-darkening was then applied to the power 
spectrum by multiplying by the inverse of R'(F) to complete the cal- 
culation. Examples of the resultant power spectra, normalized to a 
value of 10 in the interval 4 — 5 mHz for display purposes, are shown 
in Figures 6 and 7- 

The statistical accuracy of the computed spectra is estimated 
in terms of equivalent degrees of freedom and their associated con- 
fidence limits. For a one- dimensional time series for which the 
power spectrum is computed using the Blackman-Tukey (mean-lagged- 
product) method the degrees of freedom, which may be thought of as 
representing the number of estimates of power in the frequency inter- 
val Af = l/2M, are given by (Blackman and Tukey, 1958) 


where 

k = degrees of freedom, 

T = length of data set under consideration, and 
M = maximum lag in the autocorrelation function. 

A measure of the probability that an estimate falls within an upper 
and lower bound, the ratio of which is designated the confidence 
factor, is given by 
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conf 


= 10 b/(!O/k=I) = 


2.5 b 
10/k^l, 


where 

f f = confidence factor for a given confidence; 
k = degrees of freedom; and 

b = factor depending on desired confidence, given as 

follows : 


confidence 



This is an approximation, but for k ^ 4 is very close to more exact 
calculations based upon chi-square tables (Edmonds, 1966). In this 
study the least figure for k was 18.5, obtained for three of the 
records . 

It must be emphasized that no study of power spectra can be 
considered reliable or complete without the inclusion of an error 
analysis, particularly in cases where the signal- to-noise ratio is 
very poor and the peaks in the calculated power spectra are nominal. 
This must be true no matter how the power spectrum is computed, for 
the confidence limits are merely an expression of the uncertainty 
principle applied to the analysis of a time series. 

Of the 40 power spectra computed, one for each record, only 
one, shown in Figure 8, showed any evidence of having a statistically 
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significant periodicity present. The remaining spectra all showed 
structure that fell well within the confidence limits. This result 
was not unexpected, for the large antenna b'eamwidth included the 
entire disc of the sun in its pattern. Unless there were several 
exceptionally strong coherent sources, or alternately, adjacent 
regions on the surface of the sun were oscillating in phase, the 
assumption that organized activity could be confined to regions of a 

size on the order of that of supergranules meant that one was 

3 

observing on the order of 2.5 x 10 (Leighton et al. , 1962) separate 
incoherent sources, if in fact intensity oscillations did exist. 

In order to test the thesis that periodic activity at micro- 
wave frequencies is correlated with solar activity in general (Yudin, 
1968 and Simon and Shimabukuro, 1971) activity indices related to 
the amount of solar activity present in H^, soft solar x rays 

o 

(2 — 12 A), and each of several different radio frequencies were 
constructed for each data record. This was done in order to test 
for possible correlation between periodicities in the microwave 
spectrum and activity at different heights in the solar atmosphere. 
The indices were defined as follows : 




where 
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= activity index, 

T* = effective record length, and 

= apparent area of flare visible on the solar surface, 
as given in Solar Geophysical Data . 

This definition was adopted from Sawyer (1967) and is similar 
to the one used in Solar Geophysical Data since 1970 to classify 
daily H a activity. The summation in this and the following defini- 
tions includes only flares occurring during the record period. 




x-ray activity index; and 

total integrated flux above background, as measured by 
Explorer 35 (Drake et al. , 1969 ). 

I 15 = T* £ f i t i 


where 


Z x = 


F = 
x 


where 

I^t- = 15,400 MHz activity index, 

f = average integrated flux above background for 15,400 
MHz flares, and 


= flare duration . 
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The data used for determining was that taken by Sagamore 
Hill Radio Observatory and reported in Solar Geophysical Data . 

For the sake of completeness activity indices were also 
defined similar to 1^. for radio activity at 8800 MHz, 4995 MHz, and 
2695 MHz, although it was realized that there is a certain degree of 
correlation between activity at these different frequencies. There is 
also no a priori reason to expect any relation between oscillatory 
phenomena at chromospheric heights and the higher, coronal-based 
radio activity. 

Once the various indices were established for each data record 

the 40 days were divided into two groups under each type of activity 

index: those days having the highest index were defined as being 

"active" days and those with the lowest as being "quiet" days. The 

grouping under I„ , I , and for 2695 MHz had 20 days in each group, 
a x 

The grouping under 1-^ and for 8800 MHz and 4995 MHz had fewer days 
in the "active" group since flare activity was rarer for these fre- 
quencies. Consequently for these indices a record in which any kind 
of activity was manifested was defined as "active." Table 1 shows 
the list of records and to which group, by index, each was assigned. 

After classifying the records into "active" and "quiet" groups 
under each activity index the individual power spectra in each group 
were averaged to obtain an "average" power spectrum. This had the 
effect of increasing the statistical reliability of the averaged set 
over that of the individual power spectra, effectively enhancing any 


19 


periodicity that might be common to all of the spectra but masked 
because of a low signal- to-noise ratio. 

The following formula was adopted to average the spectra: 


P. 

1 



2 

L a i ex P 
d=i J 



3 


2o/iT 

j 



p. . 

1J 


where 


weigh 

weight 


P^ = i'th point in the averaged power spectrum, 

P^ = i'th point in the j ' th power spectrum, 

N = number of spectra being averaged, and 

Oj = standard deviation about the mean of data in the origi- 
nal j'th pre-whitened data set. Since the data in the 

j ' th data set were divided by a . after the mean wan 

0 

removed prior to calculation of the power spectrum, 

this weight is necessary to restore the proper relative 

units to P. . . 

ij 

kj = degrees of freedom in j'th computed power spectrum and 
b =8 for 50 ijo confidence limits. 

The exponential weighting factor is included to properly 
individual spectra according to its probable error. Since the 
of an experiment in an average is proportional to the recip- 


rocal of its uncertainty this results in emphasizing the more 
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reliable ones. The uncertainty in individual spectra was taken to 
be the probable error, or the square root of the 50 % confidence 
factor defined above. 

Confidence limits for the averaged spectra were computed by 
assuming that its equivalent degrees of freedom equaled the sum of 
the degrees of freedom in the individual spectra comprising the 


average. 
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IV. RESULTS 


A. Results of This Study 
1. Long Data Records 

The results of averaging the spectra by groups are shown in 
Figures 9 through l4 for the frequency range 1 — 15 mHz. The upper 
curve in each figure is the averaged spectrum for the "quiet" group 
under the given index, and the lower curve is for the corresponding 
"active" group. The two curves in each figure are normalized to the 
value 10 in the interval 5 — 6 mHz and separated by a decade for 
display purposes. 

It is immediately apparent that there are no peaks in any of 
the spectra that could be considered as statistically significant 
compared to the 9&la confidence limits. The fact that there are no 
outstanding peaks in either of the averaged spectra indicates that 
the antenna and receiving system do not introduce any periodicities 
in this frequency range. There does appear to be a difference in 
the general small-scale structure between the quiet and active curves 
in each case. In order to provide a more direct comparison between 
the spectra representing the active and quiet days, the ratio bet- 
ween the two spectra for a given index was taken. These results are 
plotted, by index, in the upper part of Figure 15 and provide a 
semi quantitative method of comparing the relative power in each 
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spectrum as a function of frequency. Error bars were assigned by 
adding the probable error ( 50 $ confidence) of the active and quiet 
spectra to yield the probable error for the ratio and then computing 
the 96 $ confidence limits from this result. 

An examination of these ratios shows no significant period- 
icities present in the averaged active spectra that are not present 
in the averaged quiet spectra, though several peaks are present that 
appear to be marginally significant. Peaks common to and x-ray 
activity appear at 8.0, 9-3, and 12.0 mHz and there is a general 
morphological similarity between the two curves. This similarity is 
presumably due to the fact that 13 of the records used to calculate 
the averages were classified as being active for both indices, thus 
providing a certain degree of expected correlation. A fairly well- 
defined peak common to the 8800 MHz and 4995 MHz indices is present 
at 4.3 mHz but is not considered to be statistically significant. 

It may be noticed that there appears to be a general upward 
trend, shown by light lines in Figure 15, in the ratios for the 
H Q and x-ray indices. This trend would be expected if the background 
"continuum" of the averaged active spectrum falls off with frequency 
slower than that of the averaged quiet spectrum. This result indi- 
cates that turbulant activity in the chromosphere may show up as a 
power shift toward higher frequencies in the power spectrum, even 
though no discernable periodic structure appears. The absence of a 
similar trend in the ratios for the radio frequencies is interpreted 
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as meaning that little or no correlation exists between coronal-based 
radio activity and the existence of chromospheric-based oscillatory 
phenomena . 

In order to test the possibility that the peaks evident in 
the ratios in Figure 15 were real but marginally significant, the 40 
days were subdivided into two groups of 20 days each. Each sub- 
group of 20 was then divided into 10 "active" and 10 "quiet" days 
using the activity indices as guides, and the power spectra of each 
group averaged according to the method of the preceding section. 

This averaging was done for both the I„ and I indices: the assign- 

ri x 

a 

ment of the data sets to their respective "active" and "quiet" 
categories is indicated in Table 1. The ratios of averaged active 
to averaged quiet spectra were computed as before, providing a set of 
4 curves shown in the lower-middle part of Figure 15. For each index 
a curve is labeled 68 or 69, that being approximately the year into 
which each of the 20-day groups fell. 

Inspection of these curves shows that for the two indices 
there is no significant similarity between the structure of the 1968 
averaged power spectra and that for the 1969 data. This indicates 
that the peaks in the original ratios calculated using all 40 days 
are essentially statistical in nature. It may be noticed also that 
the upward trend is in evidence in each of these ratio curves, 
especially those of H^, activity. This tends to confirm the reality 
of the trend. 
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2. Short Data Records 

An investigation was also made into the effects of calculating 
and averaging the power spectra of many short data sets rather than 
those of fewer but longer data sets. This was motivated in part by ■ 
the necessity of comparing the results of several authors (done in 
the following section), each of whom used data sets of different 
lengths. A more important reason stems from the fact that the data 
may be nonstationary in character; that is, the power spectrum of 
a data set of finite length changes with time. This nonstationarity 
would presumably result from dynamic activity in the solar atmosphere 
and the fact that few processes there remain constant for long periods 
of time. 

This is particularly true of the solar oscillations. 
Bhattacharyya (1972) in his study of velocity oscillations at several 
photospheric and chromospheric heights found that the oscillations 
exist in the form of localized bursts in which the fluctuations rise 
to a maximum value and then die down. The periods between bursts 
and the phases of the oscillations appeared to vary randomly and 
there was little oscillatory motion between bursts. It was found 
that the mean duration of the bursts decreases with height, as does 
the mean fraction of time over which the oscillations exist. These 
findings imply that velocity data, and hence the attendant intensity 
data for the upper photosphere and lower chromosphere, recorded from 
a fixed location on the sun are highly nonstationary over short data 
periods . 
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The present experiment used a large antenna beam-width, thus 
recording the fluctuations of many oscillating elements and reducing 
the signal-to-noise ratio to a near zero value. The nonstationarity 
of individual bursts would then be effectively cance lle d by the super 
position of many bursts. It is nevertheless instructive to test the 
possibility that enough nonstationarity remains such that any mar- 
ginally significant periodicities are smeared out in the averaged 
power spectra when long data records are used. This nonstationarity 
would be the result of the periodicity not maintaining a constant 
phase or frequency over the time of each individual record. 

The averaged power spectrum of several short, contiguous data 
records will not, in general, be the same as the power spectrum 
computed by considering the entire data set as a single unit. To 
see why this is so consider a data record f(t) representing the 
intensity of a fixed point on the solar disc as a function of time 
on the interval [T q , T^]. Suppose f(t) is composed of n contiguous 
data sets such that 

f(t) = f x (t) + f 2 (t) + ... + f n (t) 

where f (t) is defined on the interval [T^ T^] and is zero else- 

where . 
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f x (t) f 2 (t) f 3 (t) f n (t) 



Figure 1 

f(t) on the interval [T q , T ] (see text). 


The power spectrum of f(t) over the interval [T q , T ] is then 


PfcotM = + + ••• + ? n («>) 


n 


1 ?*(<”> 

i=l 


where f.(uj) is the Fourier transform of f.(t). 
l i 

Likewise, the averaged power spectrum composed of the average 
of the power spectra over the individual intervals is 
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The two results differ to the degree that the cross terms in 
P^ot^ 00 ) fail to cancel. These cross terms appear in the power spec- 
trum and generally express the difference between the frequency 
components of the average envelope modulating the individual data 
sets and the frequency components in the envelope modulating the 
entire data set. They may therefore be considered a measure of the 
degree of nonstationarity of the entire data set. Since we are pri- 
marily interested in the periodic structure of the individual bursts 
we may be better served by examining the averaged spectra of the 
shorter data sets; the statistical reliability will be equivalent in 
both cases. It is in this spirit that the following analysis was 
undertaken. 

Although the foregoing remarks were directed toward the prob- 
lem at hand, the same considerations should also apply to other 
problems involving power spectral analysis of nonstationary time 
series. It is sufficient for the present to demonstrate that the 
results of a power spectral analysis can be influenced by the choice 
of computing procedures, though a detailed analysis of the problem 
would involve specifying the character of the nonstationarity and 
assumptions regarding the ergodicity of the data set. We shall not 
attempt a complete analysis here. We shall instead proceed by 
assuming ergodicity as before and by using the knowledge of the 
average scale length for the bursts as reported by Ehattacharyya 
(1972) as a guide for determining the optimum data set length. 
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It was found in the above cited study that the mean burst 

duration for the high-forming H. 486l and Mg 5172 lines was lU.2 

P 

min, the line being the higher of the two forming at an altitude 

of approximately 2200 km. This is close to the altitude at which 

the 2 cm radiation is observed to originate. The mean fraction of 

time over which the H Q oscillations were observed at a fixed location 

P 

was 0.19, yielding a mean burst period of ~ 75 min. 

Since the mean burst duration is given for the H and Mg lines 

P 

considered together, and since the mean burst duration was demon- 
strated to decrease markedly with height, we may assume that for the 
Hp considered separately the duration is less than the figure cited. 

If a mean burst duration of 10 min . is adopted for the oscillations 
the mean burst period becomes ~ 50 min. 

In the present study the 40 data sets were accordingly each 
subdivided into 50 min data sets. The power spectrum was computed 
for each of these 50 min periods using the Blackman -Tukey method with 
a maximum lag of 25 min as before. Although the statistical reli- 
ability of the power spectral estimates for each 50 min data set was 
much reduced over that of the day taken as a whole, the confidence 
limits of the averaged 50 min power spectra were equivalent to those 
of the latter. Using the activity indices constructed for H Q activity 
as a guide, all of the 50 min spectra belonging to the "active" group 
were averaged using the previously defined averaging formula. The 
same was done for the "quiet" data, and the ratio of the resultant 
averaged spectra taken as before. 
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The ratio, labeled "short records", is displayed in the bottom 
of Figure 15- It is apparent that there is no statistically signifi- 
cant periodicity evident, thus confirming the results of the study 
using the full data sets. The upward trend is again in evidence, 
although not to as marked a degree as in the previous experiments 
using the H indices. The detailed structure is similar to that of 
the H ratio, confirming that the data taken as a whole is nearly 
stationary and that little or no periodic activity is present. 

We thus conclude that for the bO days of data used in this 
study there is no evidence of statistically significant periodic 
activity. Such activity may indeed exist but its amplitude must 
either be below the level of detection or it could consist of peri- 
odicities with a dispersion in periods sufficiently wide to smear out 
any discernable peaks. We also conclude that because of the absence 
of detectable constant periodic activity in this experiment there is 
no apparent correlation between solar activity and the oscillatory 
phenomena. The increased activity may nevertheless be manifested by 
the shift of power from low to high frequencies as indicated by the 
slight upward trend in the ratio curves. 

B. Comparison With the Results of Similar Studies 
Hie lack of detection of any periodic activity in this experi- 
ment was unexpected in light of similar studies made of the oscilla- 
tory phenomenon at microwave frequencies. In particular those studies 
by Yudin ( 1968 ), Durasova et al. (l97l) , and Simon and Shimabukuro 
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(1971) all reported having detected intensity fluctuations in solar 
microwave radiation. Their results, however, did not all agree with 
one another in the detected period. It is therefore worth reviewing 
their respective results with a view toward establishing the statistical 
reliability of the peaks in their power spectra and comparing, from 
a common base, those results with those of the present study. 

1. Yudin (1968) — Fluctuations at 3.3 cm 
This was the first reported observation of intensity modula- 
tion of solar microwave radiation obtained from data taken on an 
instrument with an angular resolution of 2k arc min. It was found 
that when active regions appeared on the sun radiometer output fluctu- 
ations increased markedly, with a frequency of repeatability of 
70 — 80 rfo over a period of 3 l/2 years. 

A power spectrum analysis was performed on data corresponding 
to the appearance of active regions on the sun and the results com- 
pared with power spectra obtained from data from the quiet sun. Wo 
mention is made in the paper of hew an active region is defined, but 
is taken to mean a localized enhancement of activity as evidenced by 
spectroheliograms or the appearance of plages, flares, or other 
phenomena commonly associated with solar activity. The power spectrum 
of an active period typically shewed a large peak near 300 sec with 
an intensity of 2 — 4 times that of the quiet sun (see Figure l6). 

Although it is not explicitly stated how the power spectrum 
was computed the absence of a large low-frequency (zero) peak implies 
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that the data is nearly all random noise, except for that part 
appearing in the peak at ~ 300 sec. This is incompatable with the 
assumption that intensity fluctuations result from turbulence in 
the upper photosphere and lower chromosphere. The power spectrum of 
turbulence -induced oscillations should have a large low-frequency 
component (Noyes, 1967), and since the large antenna beamwidth assures 
coverage of nearly the entire disc of the sun, including non-active 
areas, such a component is expected. This would be in agreement with 
what was found by Simon and Shimabukuro (l97l) and the present study. 
The absence of this low-frequency power suggests, therefore, some 
kind of data pre-processing, perhaps similar to that done by 
Dura s ova et al. (197I; see section 2 below) that would alter the true 
shape of the spectrum. 

Although no mention was made of how the confidence limits in 
the computed spectra were estimated one can readily calculate them 
from the information given in the paper. In this and the following 
sections it is assumed that confidence limits as calculated from 
Blackman and Tukey (1958) are the true measure of statistical reality. 

An estimate for the equivalent degrees of freedom can be made 
from a knowledge of the frequency resolution of the power spectrum 
and the length of the data set, or quantities derivable from these 
two. This is true regardless of the actual method of calculating 
the power spectrum (mean-lagged-product or fast-Fourier-transform) 
since frequency resolution and maximum frequency are set by the 
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uncertainty principle. We can therefore treat a power spectrum as 
if it were computed using the mean- lagged-product method to derive 
estimates of the statistical parameters necessary for determining 
the confidence limits. 

Figure 1 6 taken from the cited paper shows three power spectra 
of data taken from the active sun, each with a broad peak near 300 
sec. The lengths of the data sets from which the spectra were com- 
puted are given as 120 min, 70 min, and 58 min, respectively. The 
frequency resolution appears to be Af « 0.5 mHz, from which the 
equivalent maximum lag M of the autocorrelation function would be 
1000 sec. Using the formula for equivalent degrees of freedom given 
in the previous section we calculate k as being 13.7, 7-8, and 6.3 
for the three spectra. This corresponds to 96$ confidence factors of 
5.0, 9.2, and 12.4, respectively. These factors provide wide margins 
within which the respective peaks should fall in order to be statis- 
tically significant. The reported peaks all appear to be significant, 
therefore, although only marginally so. 

2. Durasova et al. (l97l) — Fluctuations at 3.3 cm 

A second measurement was made at the same frequency and by 
us ing one of the same instruments as Yudin (1968) but adding a second 
instrument for cross-correlation studies to help eliminate tropo- 
spheric influences. The effective bearawidths of the 2-antenna com- 
bination was 33 arc min, thus covering the entire disc of the sun. 
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The method of examining the data for periodic structure is 
given as being similar to that used by Yudin. We note that again 
there is no low-frequency peak in any of the spectra presented, 
which may be due to the stationarization by "... continuous summa- 
tion ..." referred to in the paper. As in the previous case, it is 
unknown how much this procedure would alter the appearance of the 
spectrum. There was no attempt by the authors to place confidence 
limits on the resultant spectra. 

The cited paper gives the maximum lag of the autocorrelation 
function as being ~ 47 min, from whence the frequency resolution 
Af ~ 0.17 mHz. An examination of the small-scale structure of Figure 
17 indicates that this must be approximately correct though the 
spectra are drawn as being continuous. The three curves of Figure 17 
are reported to have been calculated from data sets respectively 180 
min, 178 min, and 159 min long, for which the equivalent degrees of 
freedom are calculated to be 7.0, 6. 9, and 6.1. These figures trans- 
late into 96$ confidence factors of 4 . 6 , 4 . 6 , and 5«0, respectively, 
making the individual peaks in the spectra appear to be only margi- 
nally significant. Figure 18, however, shows the average of 22 
power spectra with two well-defined peaks at 280 and 400 sec, and a 
broad peak at ~ 900 sec. If we assume a figure of 7 for the equi- 
valent degrees of freedom for each individual power spectrum used in 
making the average we arrive at a 9 6$ confidence factor of 1.6. The 
peaks thus appear to be significant. 
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3. Simon and Shimabukuro (l97l)— Fluctuations 
at 3-3 and 3.5 mm 

The last of the three reports of periodic activity at microwave 
frequencies was made by analyzing data taken on the NRAO Kitt Peak 3 6 
ft radiotelescope with a beamwidth of 1.2 arc min at a wavelength of 
3.5 mm and. on the Aerospace Corporation 15 ft radiotelescope with a 
beamwidth of 2.8 arc min at a wavelength of 3-3 mm. The power spectra 
were computed using the Blackman-Tukey algorithm for the Kitt Peak 
data and the fast-Fourier-transform method on the Aerospace data. 

No confidence limits were placed on the presented spectra by the 
authors. 

One can compute the expected confidence for the various 
figures shown from data supplied in the body of the paper. Thus 
Figure 19 shows the power spectrum of a single data set of length 
50 min taken with 1.42 sec data sample spacing for a total of 2l60 
samples. Af is specified as being ~ 0.30 mHz, from which we calcu- 
late the equivalent maximum lag of the autocorrelation function as 
being M ~ 1670 sec. This is more than half the length of the data 
set so we must conclude that M is actually 1500 sec, setting 
Af = 0.33 mHz. The equivalent degrees of freedom for this power 
spectrum are then ~ 2.7 , making even the 50 rfo confidence factor on 
the order of 2. Thus the peaks shown at 3-5 and 5-6 mHz are defin- 
itely of questionable reality. 
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Figure 20 shows the averaged power spectra observed over undis- 
turbed regions for both the 3-3 and 3-5 mm data. Our analysis of the 
parameters relating sampling time and data record lengths provided 
9 6$ confidence factors of 2.2 and 2.6, respectively, for these wave- 
lengths. The 96$ confidence factors for the averaged power spectra 
observed over disturbed regions (Figure 2l) were calculated to be 
1.5 and 3.0. It therefore appears that one may safely doubt the 
reality of the observed peaks in the averaged power spectra as well 
a.s in the individual case shown in Figure 19- 

The main results of the three previous papers that are of 
interest to this study are summarized below along with an assessment, 
based on the foregoing analysis, of the reality of the principal 


peaks in the reported power 
Investigator 

spectra. 

Period 

(sec) 

Frequency 

(mHz) 

Assessed Reality 
of Results 

Yudin (1968) 

275 

3.6 

marginal 

Durasova et al. (l97l) 

4oo 

2.5 

apparently 


280 

3.6 

real 

Simon and Shimabukuro (l97l) 

180 

5.6 

doubtful 

Present work 

none 

_ _ 

real 


It is difficult to reconcile the results of these studies 
and the present one with one another. The small beamwidths of the 
instruments used by Simon and Shimabukuro should have resulted in a 
decidedly better signal-to-noise ratio with respect to the intensity 
modulations than the large beamwidths used in the other studies. If 
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the oscillations exist as coherent intensity fluctuations on a size 
the order of that of a supergranule the 3 mm instruments would be 
looking at ~ 5 — 10 separate oscillating sources, whereas for the 
centimetric data ~ 2.5 X 10^ sources would be visible. Additionally 
the millimetric study would be looking at fluctuations at a lower 
altitude than the centimetric studies. Hence, if one assumes that 
the relative amplitudes of the oscillations decrease with height, one 
would expect the 3 mm study to show a peak in the power spectrum if 
amplitude fluctuations were present. It is not immediately clear 
why large peaks therefore appear in the 3.3 cm studies but none in 
the millimetric spectra. The same might also be said concerning the 
lack of any detectable period in the present 2 cm study versus the 
large peaks in the two 3.3 cm reports; all three experiments were 
made at comparable wavelengths and all used instruments of large 


antenna beamwidth. 
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V. A MODEL FOR OSCILLATION BURSTS 


The recent work of Bhattacharyya (1972) referred to in the 
previous section suggests a model for the morphology of the oscilla- 
tions at chromospheric heights. This model may be used as a starting 
point for estimating several useful parameters that may be helpful 
in future planning for the detection of these oscillations. In the 
following sections the model is given and estimates of the minimum 
signal- to-noise ratio for which a signal can be confidently detected 
are derived. Estimates for an upper limit to the ratio of the mean 
peak-to-peak to ambient temperature for the oscillations at chromo- 
spheric heights above a supergranule are made. An upper limit for 
the mean fraction of the quiescent solar surface experiencing oscil- 
lations at any given time is also estimated, and an expression for 
the minimum amount of data needed for detection of the periodicities 
at microwave frequencies is given. 

A. The Model 

As mentioned in section IV, Bhattacharyya found the oscilla- 
tions to exist as localized bursts with well-defined times for mean 
burst duration and mean burst period. At heights characteristic of 

H_ line formation the bursts were found to have a mean duration of 
P 

14.2 min, and we set a value of ~ 50 min for the mean burst period. 
Although the physical area over which the oscillations remain 
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coherent at chromospheric heights is unknown, the work of Bhatnagar 
and Tanaka (1972) suggests that the fluctuations may be coherent 
over an area approaching that of a supergranule. 

Figure 22 shows how these bursts would appear if they are 
exponentially damped with a l/e damping time of 15 min. Although 
the period of the oscillations is shown as being 5 min, it may in 
fact be shorter. The actual period is unimportant in the following 
analysis as long as it is short compared to the damping time. Addi- 
tionally, there may be noise superimposed on the bursts. 

B. Derivation of the Minimum Signal-to-Noise Ratio for Which 
A Monochromatic Signal Can Be Confidently Detected 

In order to determine the detectability of a signal imbedded 
in noise we begin by considering the case of a monochromatic sinusoid 
of constant amplitude A impressed on stationary noise approximated 
by Gaussian noise of standard deviation o. Let 

= noise power per elementary frequency interval Af, 

B = signal power per elementary frequency interval Af, and 
s 

S = signal-to-noise ratio. > 



B 

s_ 

B m 
n 


( 1 ) 


Then 
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where m is the number of elementary frequency intervals Af in the 
power spectrum containing the full power. Now 


mAf = 


1 

2t 


where 


x 

s 


is the data sample period in seconds. 


so that 


B 

S = 2t Af 

s B 
n 


or 


B 

n 


S 

2t Af 
s 


Let 

f = 9 6$ confidence multiplicative factor 


= /( 96% confidence factor) 


= exp 


2.88 

/k - 1, 


where k is equal to the degrees of freedom. For 9 &f 0 confidence signal 


detection we must therefore have 



4o 


B + B B 

-2 1 S l + iif 

B B c 

n n 


from which 


S s 2 t Af (f - 1) 
s c 


Let denote the minimum signal-to-noise ratio for detection 


with 96 % confidence. Then 


S = 2t Af 
n s 


I 2.88 . 

exp ■ ■■■■ - 1 

l/k - 1, 


( 2 ) 


Now 


k. 2 p-J 
*>m 3 


2 — for large amounts of data (N » m), 
m 


where 

N = total number of data points in data records and 
m = equivalent maximum number of lags in the autocorrelation 


function 



4i 


= the number of elementary frequency intervals in the power 
spectrum. 

Therefore 


k 



2T 

1 

2t Af T s 
s 


= 4TAf 


(3) 


where T = Nt s , the total amount of data in seconds. Substituting 
Eq. (3) into (2) and asstiming k » 1, we obtain 


S 


n 



Af 


exp 


f— I 

/TAf , 



In the above expression for and in the following sections 96 $ 
confidence is used to define detectability. Figure 23 shows S n as 
a function of k for several different confidences with m = 50. 

Figure 24 shows the results of one of several tests where signals of 
several different amplitudes were impressed on random noise with an 
added trend. These tests confirmed the results of the derivation of 
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C. Application to the Model to Calculate the 

Burst Amplitude for the Quiescent Sun 

From Eq. (l) the minimum detectable amplitude A . of a mono- 

min 


chromatic signal is 


A . = 0/2S - 

min n 


o is the rms noise as viewed in the digitized data about the instan- 
taneous mean in a given data set. As such it depends upon the inte- 
gration time for each data sample. If a is taken as the standard 
rms noise for a data sample integration time t. of 1 sec, we have 


a 



Hence 



It is desired to calculate the minimum detectable Al/l (the 
ratio of burst peak-to-peak to ambient temperature as measured at 
the receiver output) for a single burst exponentially damped as 
shown in Figure 22. 



b3 


We define 


(power in a monochromatic signal) 

(power in a single damped burst wave train) 


1 

p [H [h(t)] 2 dt 

D JO 


where 


= burst period, 

h = waveform modulating the burst 
"G 


-t/' 


= e 


for exponential damping, and 


t, = mean burst duration, 
b 

If AI is equal to the minimum detectable peak-to-peak burst 
min 

amplitude. 


AI . = 2A . g = 2cr g 

min min o 


2S 


n 


J T i 


where it is assumed that the spreading Af^ of the peak in the power 
spectrum due to short burst duration is such that Af^ < Af . For the 
present model Af^ « 0.l6 mHz so this criterion is satisfied. 

When the sun is in the quiescent state the visible surface 

3 

contains ~ 2.5 x 10 supergranules. If each is considered to be an 



44 


independent source of coherent oscillations of random phase and mean 
peak-to-peak amplitude AI and if denotes the number of sources 
viewed simultaneously in the radiotelescope beam. 


AI . 
min 

I 

o 




where AI . /i is the minimum detectable ratio of temperature fluctua- 
min' o 

tion peak-to-peak amplitude to mean ambient temperature I q . Thus 



This equation assumes that the frequency dispersion of the bursts 

over N sources is smaller than Af. 
s 

One can evaluate these relations for the parameters of the 
present study to calculate Al ^. ^/l^. For the total "quiet" data 

N « 2.5 X 10 5 , 

s 

t =30 sec , 

s 

cr // t7 = cj 
o' 1 

= 10 K , 
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I « 10 5 K, 

o 

T « l40 h 

« 5 x 10^ sec, 

Af = 0.33 X io -5 Hz, and. 

g = 6.67 for the present model. 


Putting these figures into Eq. (4) results in a value of 
(AI /Iq) « 0.46. The large magnitude of this figure helps to 
understand why no periodicity was observed in the averaged quiet 
spectra, for no evidence exists that oscillations of such a large 
amplitude exist uniformly scattered over the solar surface. 

The same procedure may also be carried out for the best case 
for the quiet data of Simon and Shimabukuro . Thus 


N « 10 for a beamwidth of ~ 2.8 arc min, 

s ’ 

t =6 sec, 

s ’ 

0 //FT = 6 k, 

o' 1 ’ ■ 

4 

1 « 10 K, 

o ’ 

Af = 0.33 X 10 -5 Hz, 
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4 

T » 5-5 x 10 sec, and 

g = 6.67 . 

These figures yield (AI ^. ^/l^ ) « 1.5 * 10~^. 

The fact that no periodicities were observed over quiet 
regions indicates that this figure must approximately equal the upper 
limit for the mean Al/l for a supergranule . The Al/l of a localized 
area within the supergranule may of course be larger than this upper 
limit if the oscillatory motion is not spread uniformly over the 
supergranule but is instead confined to small regicns. Such a possi- 
bility is suggested by the recent work of Tanaka (1973) where it 
was found that in plages oscillating elements occupy between 2 
and 9 % of the total plage area at any given instant. The remaining 
area showed little or no periodic motion. 

Although the difference between the physical conditions of 
the chromosphere in quiescent regions of the sun and those in plages 
are such that one could not immediately assume that oscillations in 
quiescent supergranules are similarly limited in size, it may be that 
this is actually the case. High resolution Al/l measurements are not 
available for chromospheric oscillations, but investigations into 
the phenomenon as measured in the CO vibration-rotation line at 4.67 m- 
by Noyes and Hall (1972) give a Al/l of 7$ for the high photosphere. 
If this figure is assumed to hold for the higher millimetric emission 
heights observed by Simon and Shimabukuro (l97l) and the upper limit 
of 1.5 X 10 as previously calculated is taken to be the approximate 
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mean Al/l for a supergranule in a quiescent state, an upper limit of 
~ 0.02 for the mean ratio of oscillating area to total supergranule area 
is indicated. This upper limit is in reasonable agreement with the 
figure obtained by Tanaka (1973) for plage -based oscillations in V 

D. Other Useful Parameters 

From the previous section it is apparent that in order to 
study the oscillatory motions of the quiescent chromosphere at radio 
frequencies small antenna beamwidths are necessary if the total 
amount of data needed for detection is to be kept reasonably small. 

It is possible to make estimates of the minimum detectable Al/l, 
minimum beamwidth, or the minimum amount of data required for detec- 
tion as functions of relevant instrument parameters. 

For example, if it is assumed that antenna beamwidths will be 
kept smaller than ~ 20 arc min. 


N 

s 


2.5 x 10 5 X ~ 


BW‘ 


32 '/ 


where 

BW' = antenna beamwidth, arc min; and 
32 ‘ = solar diameter. 


One also recognizes that the oscillations are not, in general, mono- 
chromatic but show an associated dispersion in frequency. Suppose 
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the oscillations are slightly but uniformly dispersed in frequency 
with a dispersion width equal to pAf, where p > 1. Then from Eq. (l) 


B 

S = -g-Z— 
B n pm 


so that 


S = 2t pAf (f - l) 
n s ^ v c 

Let Rj denote the minimum detectable mean peak-to-peak source ampli- 
tude if the beamwidth. and dispersion are known. Then from Eqs. (2) 
and (4) 


R, 


4.4 


a Q gBW ' 

_ 

o 


5 PAf 

exp 

f 1.441 

- 1 

/ T i 

- 

(/TAf) 

- 


One may invert this to solve for the minimum beamwidth needed 
to detect a given R^: 


BW 


0.23 R T I 

I o 

O e 



exp 


1.44 


/TAf 





An example of how Eq. ( 5 ) might be useful is to consider an 
ideal experiment and estimate the minimum amount of data required to 
detect the oscillations. Suppose an experiment using the NRAO Kitt 
Peak 36 ft instrument at a wavelength of 1 mm is considered. We 
assume equipment and data parameters that might reasonably be expected 
to be the following: 

aperture efficiency 20 $ 

Af = 0.33 X 10 “ 5 Hz , 

t =30 sec, 

s ’ 

P = 2 , 

BW / =0.7 arc min, 

g = 6.67 , 

a J /T i - 5 K, 

R_ a 1.5 X 10"' , and 

„ h 

I D - 10 “ K. 
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The Af and T g are chosen to allow for reasonable frequency 
resolution and maximum statistical reliability in the frequency 
range 0 — 16.7 mHz. A frequency dispersion of 2Af is allowed and g 
is calculated from the present model. is assumed to be approxi- 
mately the same as the minimum detectable Al/l calculated above for 
Simon and Shimabukuro. 

Substitution of these figures into Eq. ( 5 ) yields 
T k » 4.5 X 10 5 sec = 1.25 h 

Thus only a few hours of data would be needed in this experiment to 

achieve a statistical reliability equaling that of Simon and 

Shimabukuro. We note also that with the given small beamwidth the 

number of sources N that would be viewed is ~ 1 — 2. It should 

s 

therefore be possible to accurately determine g and either R^ or the 
mean fraction of the source surface area experiencing oscillations 
at any time by examining more than the minimum required amount of 
data . 

In this chapter several useful expressions were developed 
using the statistical properties of power spectra. These expressions 
were applied to a model of the oscillation bursts for the quiescent 
sun to yield estimates for oscillation amplitude and the mean frac- 
tion of area experiencing oscillations. One could extend the analysis 
to consider active areas of the solar surface. The assumptions 
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required to do so would include the coherence length of the oscilla- 
tions and the character of the oscillation bursts. Since neither of 
these quantities are well established for active regions at chromo- 
spheric heights, the present work does not consider them. 
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VI. CONCLUSIONS 

The conclusions of this thesis are as follows: 

(1) A power spectral analysis of 285 h of 2 cm microwave intensity 
data showed no statistically significant (> 96 % confidence) 
periodicities in the frequency range 1 — 15 mHz. No correlation 
was found between 2 cm periodicities and solar activity in H^, 
x-ray, and several microwave frequencies. 

(2) A small shift of power toward higher frequencies in the power 
spectrum of the 2 cm data was found to be correlated with solar 

and x-ray activity. This power shift is attributed to in- 
creased chromospheric turbulence due to solar activity at alti- 
tudes common to H^, x-ray and 2 cm emission. 

(3) A consistent statistical analysis of previous works reporting 
evidence for the oscillations at microwave frequencies indicates 
that confidence in these previous results is marginal. 

(k) Using the statistical properties of power spectra, an expression 
for the ratio of the minimum detectable peak-to-peak to ambient 
temperature at chromospheric heights may be derived. Applied to 
a model for oscillation bursts in quiescent supergranules and 
using the most significant results to date of experiments to 
detect the microwave periodicities, this expression yields an 
upper limit of ~ 1.5- x lcf . 
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( 5 ) Based on the most significant results of microwave periodicity- 
experiments, the mean ratio of oscillating area to total quies- 
cent supergranule area at chromospheric heights is calculated to 
have an upper limit of ~ 0.02. This upper limit is in reasonable 
agreement with the results of previous direct measurements made 
of oscillating elements in H plages. 

(6) The minimum amount of solar microwave data required for detec- 
tion of the periodicities may be estimated as a function of 
equipment and data parameters. Assuming reasonable parameters, 

i 

the minimum required amount of data for an ideal experiment 
utilizing the NRAO Kitt Peak 3 6 ft radiotelescope at a wavelength 


of 1 mm is calculated to be ~ 1.25 h. 
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Table 1 


X — Active days in full study. (X) — Active days in 
cross- 68/69 study using H and x-ray indices. 


Year Day 

Mins. 

Data 

H a 

1969 251 

540 

X (X) 

252 

270 

X (X) 

255 

270 

X (X) 

256 

510 

(X) 

258 

420 


259 

360 


260 

240 


2 66 

420 

X, (X) 

269 

500 

X (x) 

271 

48o 

X (x) 

273 

360 


274 

46o 


277 

240 


278 

48o 

X (X) 

281 

520 

(X) 

328 

390 

X (X) 

337 

360 


341 

300 


342 

240 


343 

300 


344 

420 

X "(x) 

349 

390 


1969 4o 

250 

X (X) 

121 

500 


132 

600 

X (X) 

134 

48o 

X (x) 

136 

46o 

X 

139 

510 

X (X) 

142 

380 

X (X) 

144 

360 

X (x) 

148 

420 

X (X) 

150 

420 


152 

540 


154 

710 

X 

178 

630 


185 

510 


191 

600 


193 

540 

X (X) 

195 

360 

X (x) 

211 

400 



X ray 15,400 8800 4995 



XX !*! txjjxixXXIx: XX XXX XX XXX 
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Figure 1 
Figure 2 

Figure 3 

Figure 4 


Figure 5 


Figure 6 


Figure 7 
Figure 8 


LIST OF FIGURES 

f(t) on the interval [T , T ] (see text). 

o’ n 

2 cm radiometer calibration curve, antenna temperature 
versus VCO frequency. 

Block diagrams showing the 2 cm data recording and A/D 
conversion systems . 

Typical plot of antenna temperature versus time with 30 
sec time resolution. Data records used to compute power 
spectra were chosen to exclude the end pieces of each 
data day. 

Low-pass and high-pass filter functions R(f) and R^f) 
used to pre-whiten the data prior to calculation of the 
power spectra. Post -darkening was achieved by multiplying 
the resultant spectra by l/R^f). 

One of the 40 normalized power spectra computed from the 
long (> 4 h) data records. In this and the following 
figure the statistical nature of the small-scale structure 
of the power spectrum is evidenced by its presence over 
the full frequency range with an amplitude comparable to 
that of the confidence limits. 

A second power spectrum of one of the long data records. 
Power spectrum of the single record containing a statisti- 
cally significant peak near 4 mHz. No outstanding activity 
was reported in Solar Geophysical Data for this day. 
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Figure 9 

Figure 10 
Figure 11 
Figure 12 
Figure 13 
Figure l4 
Figure 15 


Figure l6 


Averaged power spectra for quiet and active days as de- 
fined using the H activity indices described in the text. 
Averaged power spectra using x-ray activity indices. 
Averaged power spectra using 15,400 MHz activity indices. 
Averaged power spectra using 8800 MHz activity indices. 
Averaged power spectra using 4995 MHz activity indices. 
Averaged power spectra using 2695 MHz activity indices. 
Ratios of averaged active to averaged quiet power spectra. 
The top six curves show the ratios, by index, of the 
curves in Figures 9 through l4. The two sets of curves 
labeled "H^" and "x-ray" each have curves labeled "68" 
and "69". Each of these latter curves represents the 
ratio of the averaged power spectra of 10 active to that 
of 10 quiet days. The bottom curve is the ratio of active 
to quiet (as defined by H a activity indices) power spec- 
tra computed from the short (50 min) records into which 
the long data sets were divided. 

Power spectra of three daya segments given as evidence for 
the existence of the 300 sec periodicities by Yudin (1968). 
The 96$ confidence limits, computed on the basis of the 
analysis presented in the text, were added to Yudin's 
original figure. The lengths of the respective data sets 
are given as 120 min (cross), 70 min (circle), and 58 min 
(square) from which the confidence limits were computed. 
(After Yudin, 1968. ) 
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Figure 17 The three power spectra reported by Durasova et al. (1971 ) 
that show detailed structure . The frequency resolution is 
calculated to be 0.17 mHz and the 9 6 $ confidence limits 
were added to the original figure. (After Durasova et al. , 


1971.) 


Figure 18 


Figure 19 


Figure 20 


Figure 21 


The average of 22 power spectra similar to the three shown 
in Figure 17. The confidence limits were added to the 
original figure. (After Durasova et al. , 1971.) 

Sample power spectrum as given in Simon and Shimabukuro 
(1971) for a single data set of 50 min length. The 50 ^ 
confidence limits for the degrees of freedom k ~ 2.7 
indicate that the spectral estimates for this single data 
set should be considered unreliable. (After Simon and 
Shimabukuro, 1971.) 

Averaged power spectra observed over undisturbed regions 
by Simon and Shimabukuro (1971 )• The 5.5 mm curve is 
the average of the power spectra of 9 records, while the 
5.3 mm curve is the average for 13 records. The confidence 
limits were added to the original figure. (After Simon 
and Shimabukuro, 1971. ) 

Averaged power spectra for disturbed regions. The 3*5 mm 
curve is the average of the power spectra of 7 records, 
while the 3.3 mm curve is the average for 39 records. The 
confidence limits were added to the original figure. 

(After Simon and Shimabukuro, 1971. ) 
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Figure 22 The model for oscillation bursts suggested by the work of 
Bhattacharyya (1972). The figure shows how the brightness 
I of a localized oscillating source would appear if the 
bursts occur as exponentially damped wave trains of random 
phase with mean burst duration (damping time) of 15 min 
and mean burst period of 50 min. 

Figure 23 The minimum signal-to-noise ratio required for detection of 
a monochromatic sinusoid as a function of degrees of free- 
dom k for several confidences. The maximum number of lags 
m in the associated autocorrelation function used to cal- 
culate these curves was 50 , the same number that was used 
to calculate the power spectra in this study. 

Figure 24 Results of one of a series of computer tests performed to 
verify the equation for minimum signal-to-noise ratio S n > 

A signal of amplitude A was added to computer-generated 
Gaussian noise with standard deviation o. A small para- 
bolic trend was also added to the data to produce a low- 
frequency peak in the power spectrum similar to the peak 
observed in the actual power spectra of the 2 cm data. 

The power spectrum was then calculated for each of several 
different signal amplitudes. It can be seen that as the 
signal-to-noise ratio S approaches the peak correspond- 
ing to the signal becomes lost in the statistical "noise" 


of the power spectrum. 
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